Tomonaga-Luttinger liquid parameters of magnetic waveguides in graphene 
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Electronic waveguides in graphene formed by counterpropagating snake states in suitable inho- 
mogeneous magnetic fields are shown to constitute a realization of a Tomonaga-Luttinger liquid. 
Due to the spatial separation of the right- and left-moving snake states, this non-Fermi liquid state 
induced by electron-electron interactions is essentially unaffected by disorder. We calculate the in- 
teraction parameters accounting for the absence of Galilei invariance in this system, and thereby 
demonstrate that non-Fermi liquid effects are significant and tunable in realistic geometries. 

PACS numbers: 71.10.Pm, 73.21.-b, 73.63.-b 



I. INTRODUCTION 

One-dimensional (ID) electron systems can nowadays 
be studied in different material systems, e.g. by deposit- 
ing negatively charged metallic gate electrodes on top of 
a 2D electron gas (2DEG) in semiconducting heterostruc- 
tures, thereby depleting the 2DEG to form the de- 
sired structure,— or in single-wall nanotubes (SWNTs).'?- 
Such ID quantum wires have been argued to realize the 
non-Fermi liquid behavior of a Tomonaga-Luttinger liq- 
uid ('TLLj^i^i^i^ii arising as a consequence of electron- 
electron (e-e) interactions. Experimental signatures of 
TLL behavior include non-universal power laws in cer- 
tain transport properties related to the tunneling density 
of states, but many other observables also may reflect 
the non-Fermi liquid properties of a TLL. Experimental 
observations in semiconductor quantum wires^ were ex- 
plained by TLL parameters of the order of gc « 0.4 to 
0.5, while the corresponding parameter in SWNTs^ was 
reported as gc+ ~ 0.16 to 0.3. Both values are signifi- 
cantly smaller than the respective noninteracting value, 

5 = 1- 

Very recently, graphene monolayersiSili have become 

available as a new realization of a 2DEG, albeit with 
properties strikingly different from their semiconducting 
counterparts. The kinetic energy of graphene close to 
one of the Dirac {K, K') points is described by a two- 
component chiral Dirac- Weyl Hamiltonia n^^d'^ 



H — vpT • (p — A) 



(1) 



of massless relativistic particles moving at graphene's 
Fermi velocity v-p w 10^ m/sec, instead of the usual 
Schrodinger Hamiltonian p'^/2m* (with effective mass 
m*). In Eq. ([T]), r denotes the vector of Pauli matrices in 
sublattice ("pseudo-spin") space, while the physical spin 
as well as the valley {K — K') degrees of freedom are left 
implicit. Furthermore, we have allowed for a static in- 
homogeneous orbital magnetic field perpendicular to the 
graphene plane (the field components in the plane do not 
affect orbital motion), B = B(x,y)ez, which is incorpo- 
rated by minimal (Peierls) coupling in terms of the cor- 



responding vector potential A{x,y). This gives rise to 
interesting magnetic barrier and magnetic confinement 
effects^^ii^^ 

It is well known, both theoretically^^ and 
experimentallyfi^ that a magnetic-field gradient can 
give rise to unidirectional ID snake states. Such orbits 
were recently studied theoretically in graphene22i2ii22 
and in SWNTs.^- Snake states carry current along the 
lines where the magnetic field changes sign, and hence 
is zero. On a classical level, they can be understood 
as half-orbits of different circulation sense (for B > 
and B < 0), patched together to form a unidirectionally 
propagating orbit.— Pairs of snake states running 
antiparallel to each other are referred to as double 
snake statesi^^ In many regards, double snake states 
correspond to the standard left- and right-movers in 
ID quantum wires. For example, they should exhibit 
quantized conductance in multiples of 4e^//i (including 
spin and valley degeneracies). Since the snake states are 
spatially separated, this quantization should be robust: 
shallow impurities are not expected to cause scattering 
between snake modes of opposite directionality. 

In this work, we address consequences of the long- 
ranged but ultimately screened e-e interactions within 
and between the counterpropagating snake orbits in 
graphene magnetic waveguides. For a wide class of exper- 
imentally relevant field profiles, we show that a TLL state 
with broken Galilei invariance and extremely weak dis- 
order sensitivity can be realized. On a general level, the 
importance of e-e interactions for the correct interpreta- 
tion of experimental data in graphene has recently been 
stressed:^ Theories describing e-e interaction effects on 
the transport properties of electrons in graphene have 
been proposed for strong homogeneous magnetic fields^! 
and for zero magnetic field. ^S A recent debate has dis- 
cussed the question whether interacting electrons in un- 
doped graphene form a Fermi liquid or not i^^'^" Interac- 
tions are also predicted to yield a TLL state in special 
graphene nanoribbons with armchair edgesi^ 

For the related case of interacting metallic SWNTs, the 
effective low-energy theory predicts a four-channel TLL 



statCfSIiM where the spin and valley degrees of freedom 
give rise to the four channels. There is one charged (c+) 
channel, where the long-ranged e-e interactions play a 
crucial role, while the three neutral channels are basi- 
cally insensitive to interactions. We will see that the 
situation in a graphcnc magnetic waveguide is similar, 
and the parameter g discussed below plays the role of 
the SWNT parameter gc+- Albeit the interaction does 
not spoil conductance quantization in dc transport for 
adiabatically connected reservoirs;^ it nevertheless de- 
stroys the Fermi-liquid character of the system. In fact, 
it leads to non-universal power laws in the tunneling 
density of states, and to peculiar ac transport and shot 
noise^ properties at low temperatures. These phenom- 
ena are appropriately described by TLL theory.— i^i^i^ 
The respective power-law exponents can be directly in- 
ferred from Refs. |32||33| by simply replacing gc+ with our 
estimate for g, see Eq. (15^ below. 

After introducing the model and the magnetic field 
profiles in Sec. Ill A[ the bandstructure is studied in Gaus- 
sian approximation in what follows in Sec. Ill Bl The ana- 
lytical bandstructure results are validated by comparing 
to exact diagonalization results. The numerical diago- 
nalization is briefly discussed in the Appendix. The lin- 
earized bandstructure for a double-snake state waveguide 
leading to TLL behavior is then described in Sec. IIII Al 
The physics of a TLL is governed by a dimensionless in- 
teraction parameter g, for which general expressions in 
terms of certain velocities are derived in Sec. IIII I3l These 
velocities are obtained from perturbative expressions for 
the ground-state energy, and yield the analytical results 
for g given in Sec. lIVI We conclude in Sec.|Vl In Secs. lIIII 
and IIV| to be specific, we focus on electron-like excita- 
tions by assuming a positive value of the Fermi energy 
Ef- Below, we often take units such that h = vp = 1. 



II. MODEL AND BANDSTRUCTURE 



A. Model 

In this paper, we consider magnetic fields B — B{x), 
guiding particles homogeneously along the j/-direction. 
This implies that the wavenumber k along this direction 
is conserved. Two-component eigenstates of the Dirac- 
Weyl Hamiltonian llj can then be written as ip{x,y) ^^ 
G^^^i4>kix),Xk{x))'^ ■ The vector potential can be chosen 
as A = A(x)ey, with the spinor components obeying 



nent (j)k, 



-iO, 





-l-iA;- 



iA{x) 



-idx 



-\k 


Efc 



- iA(x) 

(t>k{x) 
Xk{x) 



4>k{x) 
Xk{x) 



(2) 



Complex phases may be chosen such that 0^ is real and 
Xk purely imaginary. After squaring, Eq. (U) can be cast 
into a Schrodinger-like form for the upper Dirac compo- 



[-dl + [fc - Aix)f - B{x) - el] M^) = . (3) 

A similar equation holds for the lower component Xk{x), 
with the sign of the "pseudo-Zeeman" term ~ B reversed. 
Unless Efc = 0, Eq. ^ imphes 



(4) 



dx\M^)\ = / dx\xk{x)\' = l/2. 



Note that in path-integral approaches to relativistic 
quantum mechanics, in order to guarantee convergence 
of the Wiener measure;^ often the square of the Dirac 
Hamiltonian ([T]) is considered. Path-integral represen- 
tations, on the other hand, allow for systematic ap- 
proximations, and therefore Eq. (I3|) is a useful starting 
point for the Gaussian approximation, cf. Sec. Ill Bl where 
the boundedness of the differential operator appearing 
in Eq. ^ is exploited for either sign of B{x), in the 
spirit of a saddle-point approximation. Our method dif- 
fers from WKB-type approaches recently put forward to 
describe the electronic properties of graphenei^^*^ We 
note in passing that massive Schrodinger particles obey 
a related equation as Eq. ([3]), with quadratic momenta 
multiplied by 1/(2to*), in the same magnetic field pro- 
file; only the pseudo-Zeeman term must be removed and, 
of course, the energy e^ must be replaced by et, i.e. hole 
and zero-energy states both disappear. As a consequence, 
most of our conclusions also apply (at least qualitatively) 
to magnetic waveguides based on traditional Schrodinger 
fermions. 

We consider the class of magnetic field profiles given 
by 



B{x) = lyaJBiy/tOBxy 
In our gauge, we thus have 



Bn 



A{x) 



9 V 

UJ-a X 



BnX 



(5) 



(6) 



The index v can describe rather different situations, but 
we are only interested in v being a natural number. For 
instance, for i/ = 1, we recover the homogeneous mag- 
netic field case, giving rise to the standard relativistic 
Landau levels. For i^ = 2, the profile ([5]) instead de- 
scribes a setup with one snake state propagating along 
the y-direction, while i' — S (or, more generally, all odd 
u > 1) can give rise to a double snake state geometry, 
where the background magnetic field —Bq allows for lines 
with _B = 0, and ujb sets the inhomogeneity scale. Equa- 
tion Q manifests the electron-hole symmetry e^ <-> —Sk 
of Eq. ([2]) , with a zero-energy eigenstate {sk = for all 
k) appearing whenever ly is odd, but not for even i/.^^ 

Equation ^ qualitatively describes many situations 
of experimental relevance, where typically smooth mag- 
netic field profiles are present. Of course, far away from 
the waveguide defined by the snake states, the actual 
profile is different in practice, but this does not signifi- 
cantly affect the TLL discussed below. In fact, we have 



also analyzed step-like field profiles, such as the ones de- 
scribed in Ref. 21, with very similar results and conclu- 
sions. For v = 3, counterpropagating snake states are 
centered around x — ±d with B(±d) = 0, leading to 



d^ 



^ 



V3 



LJB 



(7) 



such that 2d is the parallel distance between counter- 
propagating snake states. In this configuration, a TLL 
can be realized, and most of our analysis will deal with 
this case. 

Below, we will ignore the Zeeman splitting due to the 
interaction of the true electronic spin with the magnetic 
field creating the waveguide. A simple estimate for a 
homogeneous magnetic field already shows that this ap- 
proximation is justified in graphene. The physical Zee- 
man splitting Az = gcfisB amounts to 0.116 meV for 
B = 1 Tesla, taking g^ = 2 and the free electron mass rrio 
going into Bohr's magneton /is- This value can be com- 
pared to the orbital splitting Aorb between subsequent 
levels — in the language of Eq. ([3]), this corresponds to 
a "pseudo- Zeeman splitting" — , with the result 



(8) 



predicting that even for ep — I eV, the Zeeman splitting 
is 50 times smaller than the orbital splitting. In view 
of the smallness of the ratio ([5]) , we will neglect Zeeman 
terms in what follows. In any case, their effects on the 
low-energy theory of interacting electrons in graphene 
waveguides are standard, and could be included along 
the lines of Refs. \MMU- 



potential appearing in Eq. ([3]), 

which depends on the sublattice component r^ = ±1 
of the wavefunction. Obviously, for any k, V^^2 is in- 
variant under the combined operation r^ —> —Tz and 
e ^ -e, so that \M0\ = \Xk{-0\ in Eq. ©. For 
K — > — oo, the minima of V5/=2 approach k^ at ^ = 0, so 

that |0fc(O| = \XkiO\ ~ e-^^^«' and £fe^_oo ^ -k. 
This indicates that the velocity reaches, up to sublead- 
ing corrections, the negative of the Fermi velocity. This 
is precisely the snake state, with both components of 
the Dirac spinor localized near x = 0. In the other 
limit, K —f +00, the minima of Vi,=2 approach —2^^ at 
C = Tz\/k. In that case, the two sublattice components 

\4>kiO\ = \Xk{-S,)\ "" e"^K~^^ are spatially separated 
from one another and from the snake mode, provided the 
distance exceeds the widths of these distributions. This 
result suggests an interesting application as a "sublat- 
tice filter", where the magnetic field leads to a spatial 
separation of particles located on different sublattices. 
However, for the other K point, the sublattice states are 
exchanged, and in order to see such an effect, one would 
need to have a valley-polarized system (i.e. a single K 
point). The corresponding energy in Gaussian approxi- 
mation is Efc^+oo —^ 0. This result cannot be recovered 
using WKB-type approaches22iS which are more suited 
to describe higher excited states. Interestingly, the posi- 
tions ^ — ±y/K of these states remain protected against 
the pseudo-Zeeman field "inclination" from Eq. ([3]) , con- 
trary to naive expectation and in contrast to the non-zero 
shift found in any of the excited states. Finally, also ex- 
cited energies can be estimated in this way. For example, 
the first excited level is expected at Sk — 2uj-q k^/^. 



B. Gaussian approximation 



Next we describe our analytical approach to the band- 
structure and the eigenfunctions. They allow for closed 
form expressions of the TLL parameter g in Sec. IIVI At 
large |A:|, when anharmonic contributions of the effec- 
tive potential appearing in Eq. ^ are suppressed, the 
Gaussian approximation becomes exact. To confirm the 
accuracy of the analytical results, we have carried out 
numerical diagonalizations of the matrix representing the 
Schrodinger-like Hamiltonian ([3]) in a complete basis set. 
This is briefly described in the Appendix. 

It is instructive to first study one of the two coun- 
terpropagating snake modes individually. We therefore 
set Bq = and J/ = 2 in Eq. dS]), i.e. B{i) = 2tJBC 
with uj-Q > 0, where we introduce dimcnsionlcss lengths, 
f — ^lob, X, and momenta, n = kj ^pJjQ. The Schrodinger 
version of this model has been studied previouslyii^ The 
single snake state is now centered near a: = with (posi- 
tively or negatively charged) particles running in the neg- 
ative j/-direction. The spectrum in this case is not sym- 
metric, Efc 7^ e_fc. We then need to discuss the effective 



In effect, we then arrive at a picture where snake (near 
X = 0) and "bulk" modes (near x = iojg vfc) will de- 
velop. We here distinguish "snake" and "bulk" modes by 
their respective group velocities. Fig. [T] clearly demon- 
strates how the eigenstates evolve from snake states (at 
k — > — cx) with QfeEfe = —1) into bulk modes at fc — > -t-oo. 
In Fig. [5] the corresponding metamorphosis is displayed 
of a (single) snake state at sufficiently negative k centered 
around ^ = 0, see Fig.[21[a), into a bulk state, residing in- 
creasingly far away from ^ = with increasing fc > 0. For 
K = 8, cf. Fig. Hfb), the state is located near ^ — 2.83r2 
as expected. 

Let us now turn to the magnetic field profile with v = "?> 
in Eq. (O, which should exhibit double snake states of 
opposite directionality due to the existence of two zeros 
of B{x). To some extent, we now have two copies of the 
above single-snake state situation at a distance 2d, cf. 
Eq. ([7]). Since B{x) = B{—x), the dispersion relation 





FIG. 1: (Color online) Electron-like energy eigenvalues Ek ver- 
sus momentum k (both in units of y/tUB), as obtained by 
numerical diagonalization of Eq. ([3]) for a magnetic field pro- 
file with 1/ = 2 and Bo = in Eq. (O. The solid (black) 
curve denotes the lowest positive-energy eigenstate, and the 
dashed (red) and dash-dotted (blue) curves give the next two 
pairs of excited states. The dotted curve indicates the limit- 
ing snake-state dispersion Ek = —k for A; < 0, and the result 
in Gaussian approximation for the first excited energy band, 
Ek ~ 2LJg fc^'^, for fc > 0. Note that for the lowest state, 
Ek^+oc approaches zero energy,—!^ in agreement with the 
Gaussian approximation. 



is now symmetric, £k = £-k- At fc — > ±cx3, we thus 
anticipate the coexistence of snake and bulk modes. In 
addition, an exact zero-energy eigenstate, Sk = 0, must 
now occur as a consequence of the index theorem;^ The 
potential entering Eq. ^ is 



v,=3{K,o^ie'bo^~^f 



,(3^2 -6o), (9) 



where bo — Bq/cob- This potential is invariant under a 
simultaneous sign change of k and ^, transforming left- 
into right- movers and vice versa. Similar as for i/ = 2, 
we can obtain energies in Gaussian approximation. The 
lowest energy, only for r^ = -f 1, equals zero at large |A:|. 
This describes the zero-energy state, present for any k at 
odd i^, cf. Sec. Ill Al and the Appendix. The first excited 
(positive) energy is approximated as 



£k 



V6|wBfci'/'+5o/[V6|wBfcr/'] 



(10) 



This result is included to the numerically obtained spec- 
tra, see Fig. [3] The corresponding eigenstates are local- 
ized around^ = (JkI^/^ -h 6o/3|k|^/^ -f r^/3|K|)sgn(K), i.e. 
increasingly deep in the system's bulk with increasing 
|fc|. Fig. [3] shows a typical spectrum obtained from nu- 
merical diagonalization. Depending on the slopes dk£k at 
large |A:|, one type of bands goes like Sk — ±WF|fc|, corre- 
sponding to the counterpropagating snake states centered 
around x = id, see Eq. ([7]). Snake states move with the 
Fermi velocity of graphene at large |fc|, irrespective of the 
magnetic field profile. ^^^s^S The other bands in Fig. [3] 
exhibit smaller slopes at large |A:|. The corresponding 
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FIG. 2: (Color online) Numerical diagonalization results for 
the probability densities of the spinor components |';/>k(^)| 
(black solid curve) and |Xk(OP (^^d dashed curve) of the 
lowest eigenstate to positive energy of Eq. ^ with u = 2 and 
Bo = 0. (a) is for momentum k = —2, and (b) for k = 8. 



states are localized increasingly further away from the 
center of the wire at increasing |fc|, and we thus call them 
again "bulk" modes. Due to their different slopes, bands 
of different types should cross, and Fig. [3] indeed reveals 
avoided intersections. Such avoided level crossings can be 
attributed to some residual hybridization between snake 
and bulk modes. They become successively less impor- 
tant as the bulk state's center moves away from the snake 
state with increasing |fc|. 

Fig. Invalidates the behavior of the two types of eigen- 
states for the double-snake i/ = 3 situation with k = — 6. 
In view of Fig. [31 this value of \k\ is just beyond one 
of the (narrower) avoided crossings, so that the lowest 
positive-energy state, exhibiting a moderate slope, must 
be identified as a bulk state. This is indeed confirmed 
in Fig. m where the densities of both components of this 
state are seen. Gaussian approximation predicts its cen- 
ter to be at ^ = —2.1 + O.OBt^, which is nicely confirmed 
by comparing to exact diagonalization results. The first 
excited state at k = —6 has dkSk — —1, and therefore is 




FIG. 3: (Color online) Same as Fig.[l]but for z/ = 3 and Bo = 
1.526lob- Full curves (in different colors) give the numerical 
diagonalization results for the eigenenergies. When the Fermi 
level intersects only the lowest (black) curve, one has precisely 
one (spin- and valley-degenerate) left- and right-moving state 
in the waveguide. This leads to a TLL state. The dotted curve 
denotes Sk = |fc|, and the dashed curve gives the estimate 
(|10|) for the lowest dispersing energy band in the Gaussian 
approximation. 




|(^(5)| , lowest state 
\X&\ , lowest state 
|(^(^)| , first excited state 
|X(^)I , first excited state 










FIG. 4: (Color online) Same as Fig.[2]but for a magnetic field 
([5]) with v = 3 and Bo ~ 1.525ujb- Shown are the two lowest 
positive-energy eigenstates with n — —6. 



classified as snake state. Indeed, as seen in Fig. 01 both 
pseudo-spin components of this state reside near ^ = 0.7, 
where B{^) = 0. Note that, according to Eq. ([22]) . densi- 
ties for odd v are mirror-symmetric under simultaneous 
sign reversal of k and £,, so that the k ~ +6 state is cen- 
tered at ^ = —0.7. Finally, we note that all densities in 
Figs. [2] and 0] are well approximated by superpositions 
of suitable Gaussians, in accordance with our Gaussian 
description. 



III. INTERACTION EFFECTS 
A. Waveguide model 

When considering a ID graphene magnetic waveg- 
uide with counterpropagating snake states at low energy 
scales, see Sec. IIIBl the question arises whether inter- 
acting Dirac fermions in such a waveguide belong to the 
TLL universality class^i^i^i^ As always happens in ID, 
Landau quasiparticles will be destroyed by any non-zero 
e-e interaction strength,'^& but whether the resulting non- 
Fermi liquid is a TLL remains to be shown. Indeed, this 
expectation is corroborated by the analogous situation in 
a quantum Hall bar,"*^'^^ where, as a consequence of the 
long-ranged Coulomb interaction between different edge 
states, a TLL with spatially separated left- and right- 
moving edge states emerges.*'^ 

We start with a discussion of the relevant single- 
particle bandstructure. For a magnetic waveguide with 
ly = 3 and iJo > in Eq. ([S]), the lowest positive-energy 

subband e^ ' = e_l has left- and right-going snake states 
near x — ±d, see Eq. (O, that essentially move at veloc- 
ity ztvp. We assume that the Fermi level intersects these 
states at ifep, and that all other states are energetically 
sufficiently far away. This bandstructure was discussed in 
detail in Sec. IIIBl see Fig. [31 Using second quantization, 
the kinetic energy is then described by 



i/o = V£l')clc 



k 



k^k ! 



(11) 



where c\, (c^) creates (annihilates) a Dirac quasiparticle 
with momentum fc, and spin and valley indices are kept 
implicit. The electron field operator for waveguide length 
L along the 2/-direction is thus written as 



^{x,y) 



J_ Y^ iky 



4>k{x) 

Xk{x) 



Ck 



(12) 



The crucial parameters characterizing the TLL are cer- 
tain velocities;^ In the noninteracting case, Dirac parti- 
cles move at velocity2i 



{Ty)k == 2t;Flm / dx (l)l{x)xk{x) = dkSk = Vk 



vf{ 



given by the slope of the energy dispersion, just as for 
Schrodinger particles. Assuming that ep is sufficiently 
far away from both the band bottom of e^ and from 

the next-higher energy band, we linearize e^ about the 
two Fermi points k = zt/cp, yielding velocities ±Wfep . This 
also separates right (fc > 0) from left (fc < 0) movers in 
Eq. pT|) . and implicitly defines the standard bandwidth 
cutoff around the Fermi level used in TLL theory. 



B. Tomonaga-Luttinger liquid 

Next we incorporate e-e interactions within an effec- 
tive low-energy theory. We consider the pair interaction 



potential 



1 44 



M^(X1,X2) 



KO 



IX1-X2I ^(Xi-X2)2+4Z?2 



(13) 

between electrons at coordinates Xi = (xi, j/i) and X2 = 
{x2,y2)- This form specifically accounts for screening 
by metal gates positioned at some distance D from the 
waveguide. Its strength is governed by the dimensionless 
"fine structure constant" a = e^ /[kqHv-p], which basically 
depends only on the dielectric constant kq. For typical 
substrate materials, one has values kq « 1.4 to 4.7, re- 
sulting in a « 0.6 to 2i^^i^i^ For graphene, both the 
kinetic as well as the Coulomb energy scale (~ ^/n for 
particle density n) in the same way4^ The resulting weak 
tunability of the e-e interaction strength in graphene is in 
stark contrast to the situation in semiconductors, where 
n allows to alter the relative strength of Coulomb inter- 
actions over orders of magnitude. 

When constructing a low-energy theory for interact- 
ing Dirac fermions in the double snake state waveguide 
of Sec. IIII Al the resulting ID e-e interaction processes 
can be classified as forward-scattering and backscatter- 
ing processes^i^i^ii The spatial separation of the unidi- 
rectional snake states here implies a strong suppression 
of e-e backscattering processes, where the relevant cou- 
plings are exponentially small in the parameter kpd 3> 1. 
In the following, we discuss the regime 



kpD > kpd > 1 



(14) 



where backscattering processes are negligible. This situa- 
tion is reminiscent of the SWNT case?^ where one, how- 
ever, finds only an algebraic suppression of the backscat- 
tering couplings with increasing SWNT radius. We then 
only need to include forward-scattering processes, see 
also Ref. |3^, and arrive at a four-channel TLL model. 
The three neutral sectors involve spin and valley degrees 
of freedom, and are decoupled from each other and from 
the charge sector (spin-charge separation). For not too 
strong interactions, as expected in graphene, the neutral 
sectors will remain basically unaffected by interactions, 
with their velocity parameters (almost) equal to w/cp, see 
also Ref. |43. This implies that the TLL parameters for 
the three neutral sectors are just given by the noninter- 
acting value {g = 1). In the following, we then focus on 
the charge (c-f ) sector only. 

The resulting TLL is most conveniently described by 
Abelian bosonizationj^iSiiii?^ For the c+ sector, the re- 
sulting Hamiltonian is 



H,+ = lldy [vjidySiy)]' 



-v^idyHy)? 



(15) 



with bosonic fields subject to the algebra 
[<^{y),dye{y)] = i6{y - y'). Equation ^ reflects 
the fact that density waves (in contrast to quasiparti- 
cles) are undamped in a TLL.^" With Eq. p5|) and the 
usual bosonized form of the Fermi operators,^ almost 



any observable of physical interest can be determined 
exactly at low energies and long wavelengths, where the 
TLL model applies. In general, wj < wn can deviate 
from Wfcp as a result of the repulsive e-e interactions. 
Interaction physics is thus encoded in wj and wn- These 
velocities determine the dimensionless TLL interaction 
parameter g = gc+ and the plasmon velocity v according 



toi 



V^j/UN 



\/vyv^ 



(16) 



In principle, both parameters are experimentally ac- 
cessible, e.g. through the tunneling density of stateS)^ 
momentum-resolved tunneling,— or via plasmon propa- 
gation times.— 

The velocities wj and wn can be extracted in an elegant 
and exact manner from thermodynamic relations;^^ 



vn 



AL dkl 



v.} 



4L dS^ 



(17) 



provided the fully interacting ground-state energy den- 
sity Eq/L for fixed left (— /cp M and right (fep ) Fermi 

momenta is known, where S = ikp — kp ] /2 and 

ko = (fc^+'' + fc^"' j /2. The derivatives in Eq. dTT]) are 

evaluated at 6 = and kp = k-p, and the energy Ep in- 
cludes the spin and valley degrees of freedom. Clearly, v^ 
is proportional to the compressibility, and when Galilei 
invariance is realized, wj = Vk^ is unchanged by interac- 
tions. However, as we shall see below, this symmetry is 
not obeyed here due to the periodic superstructure im- 
posed by the snake orbit. 



IV. TLL PARAMETER 

Unfortunately, exact results for Eq are known only for 
a limited number of integrable models, such as the Hub- 
bard model^" or the Sutherland model. ^^ Even then, one 
still has to numerically solve coupled pairs of integral 
equations to access E'oj a-nd hence the velocities f j_n in 
Eq. (J17p . In actual calculations of wj and wn for noninte- 
grable models (which is the case here), one has to resort 
to approximations4^i^ 



A. Perturbation theory 

We now use perturbation theory to obtain Ep , and thus 
the velocities P7|) . for relatively weak interactions. This 
approximation is almost exclusively used in the literature 
in order to obtain estimates for the TLL parameters of 
generic interacting ID fermion systems. In that case, the 
ground-state energy can be split into three terms. 



Eq = £^kin + -^Ha 



-E^Fock , 



(18) 



from which several contributions to the susceptibilities 
(fT7)) follow. Accounting for spin and valley degeneracy, 
we have 



52^k: 



(19) 
as expected for noninteracting quasiparticles with veloc- 
ities ±Ufep. 

From the Hartree interaction term, we then obtain 



d^E^, 



L 



dkl 



r]{kF,kF)+r]{kF,-kF)+n , (20) 



where, including again both spin and valley components, 
we define 

rj{k,k') = ri{k',k) = 8 dx dx' nk{x)nk'{x') 



X aln^l + ^ID/ix'-x')]^ , (21) 

with the particle density for wavevector k at location x, 

nk{x) = \<pk{x)\^ + \Xk{x)\^ - n.ki-x) . (22) 

Note that J dx nk{x) = 1, see Eq. (|3]). Furthermore, we 
have introduced the quantity 



n 



fcp 



dk dkp [r]{k, kp) + r]{k, -kp)] 



= 2/ dkdkp[ri{k,kF)+ri{k,-kF)] 
/o 



(23) 



The second equalities in Eqs. (|22|) and p3|) are valid for 
all odd V > \. In Eq. (HI]), we have used that only long 
wavelengths q — *■ are important in i^Hartrcc- To see this, 
consider the one-sided Fourier transform of Eq. p^ , 



W{q, x) = a dy e'' 



,igv 



1 



1 



a 



(24) 



K„{\qx\) - Ko{\q\V^D^+^) 



implying that W{q -> 0, a;) ~ a \n y/T+\2D/x^ . With 
these definitions, we obtain in a similar manner 



d i^Harticc 



L 



ri{kF,kF)-ri{kF,~kF)+n . (25) 



In view of Eqs. (fT6|) to (f25|) . only the magnitude of 
r]{kF,—kF) yields a nontrivial contribution to g at this 
level of approximation. Within the g-ology terminologyj^, 
we may identify this term with 32, measuring the strength 
of forward scattering between particles of different direc- 
tionality (henceforth referred to as chirality). On the 
other hand, scattering between equal-chirality particles 
is described by g^, identified here as r]{kF,kF)- How- 
ever, according to Eqs. PP)) and PSJ) . an additional 
term TZ is present, which effectively modifies (74. In 



ID quantum wires with continuous (Galileian) transla- 
tional invariance, both n^^x) and ri{k,k') are indepen- 
dent of k.k' , and hence TZ = and 3|£'Haitioc = 0. In 
that case, the TLL parameter g depends only on the 
zero-momentum Fourier component of the interaction, 
W{q ~ 0) w a\n{D/d), where d is of the order of the 
wire width. 

Similar (though slightly more involved) expressions can 
be found for the Fock contributions to Eq. (|17p , which are 
of the order W{q ~ 2fcF). Using Eq. p^ . this amphtude 
can be estimated as 



W{2kF 



^—ikpd 



Skpd 



which is parametrically smaller than the Hartree am- 
plitude. Similar to backscattering contributions, Fock 
contributions can thus safely be neglected against the 
Hartree terms for the parameter regime (J14p . 



B. TLL parameter estimate 

In order to estimate the magnitude of the above terms, 
in particular of ?y(fcF, — ^f), it is necessary to have some 
handle on the unperturbed wavefunctions (pkix) and 
Xk{x), together with the resulting densities nk{x). We 
approximate their density profiles as 



n.iO 



{l2boK^ 



'(126oK^)^^''K+sgn(K)v'JJ5-c 



(26) 



Note that the true densities describing snake orbits, see 
Fig. m are somewhat more complicated, with a double- 
peak shape. However, the simplified single-Gaussian 
form in Eq. (|26p captures the essential physics and al- 
lows for analytical progress. 

Given Eq. ([26]) . introducing the lengthscale 



A 



-Bow|fcj 



(27) 



and using Eq. (j2ip . we are now in a position to estimate 



r]{kF,-kF) 
r]{kF,kF) 



Baln(L>/d) , 



an 



(e'^'^D/x) 



(28) 



assuming I? > d > A, see Eq. pi)). Here C = 
0.577 ... is the Euler constant. We see that the ratio 
ri{kF, —kF)/ri{kF, kp) approaches unity for D ^- 00 as in 
Galilei-invariant ID wires, where in addition one also has 
TZ = 0. In order to estimate TZ, we first observe that in 
Eq. (P5)) . the 77(fc, — fcp) term is suppressed by a factor 
g-8(d/A) g^g compared to ri{k,kF). This factor becomes 
small in the regime (|14p . and we can therefore approxi- 
mate TZ c^ 2 L dk 9fep77(fc, /cp), which can be calculated 
in closed form for the density profile (pS)) . 



TZ 



8a 



V2(ci + 2c2) + 4(%/2ci - 1) ln(i:>/A) 



a[5.73-2.601n(i:)/A)] 



(29) 



Employing the incomplete Euler Beta function, we find 
ci = i UV2 -Jdt ^^^) - 0.6496 (30) 



and 



rl ^1/4 

C2 = / dt In 

Jo VTTt 



1 + t 



l + Vi 



-0.07171. (31) 



Remarkably, TZ decreases with increasing D/X, changes 
sign at D ~ 9. 02 A, and then continues to decrease loga- 
rithmically. Although asymptotically smaller by a factor 
4(^/2cl - 1) ~ 0.326 than the leading contribution ^ to 
34, the 7?.-term is important for quantitative estimates of 
the TLL parameter. For example, when 7?. < 0, standard 
expressions (without TZ) would overestimate g, pretend- 
ing too weak interaction effects. The usual expressions 
are based on Galilei invariance, which is broken in the 
present system due to the periodic superstructure im- 
posed on the ID wire by the snake orbits. It is not ob- 
vious to us how standard estimates to 174, starting from 
the microscopic interaction P^ . would recover the 7?.- 
contribution. 

Combining Eqs. (fTB)) to (P^ . we get the analytical es- 
timate for the TLL parameter. 



9 



go - InjD/d) 
go + HD/d) 



1/2 



(32) 



in the regime D ^ d ^ A, see Eq. P^ . where A is given 
in Eq. (p7| . Here, we have abbreviated 



9a 



2a 
In 



\/2(ci + 2c2) + 4(\/2ci - 1) In 



D 



e'^/2 — 
A 



(33) 



The TLL parameter (p2)) is depicted in Fig. [5l where 
It^feSfepI = ''^F and a "fine structure constant" a = 1 have 
been assumed. However, according to Eq. ([55)) . changes 
in a may be compensated for via changes in A, i.e. by 
modification of kp or of the magnetic field parameters. 

As shown in the inset of Fig. [3 values of about g ^ 0.3 
are expected for the chosen value of Bq, with no pro- 
nounced variations when changing D/X. Such weak sen- 
sitivity to kp found in graphene is in stark contrast 
to semiconducting wires, where g can vary significantly 
when changing the carrier density4^ As discussed in the 
Introduction, similar values for g as compared to the val- 
ues in Fig. [5] have been reported for other TLL systems 
such as quantum wires^ or SWNTsi^ '^^i'^^ On the other 
hand, the main panel in Fig. [5] demonstrates that g can 
be widely tuned in graphene wires by changing the snake- 
state separation 2d in Eq. ([7]) via the magnetic field pa- 
rameters, in particular by sweeping the background field 
Bo. 
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FIG. 5: TLL parameter g in Eq. (|32j for a = 1 and D = 

IOO/y^ojb. Main panel: g as a function of Bq/uj-b, for fcp = 
IQy/iUB. Inset: g vs kp/^/ZuB for Bo = 1.526ii;B. Note that 



the regime X <^ d translates into kp S> -^lubBq 
here implies k-p/^/coB > 1. 
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V. CONCLUSIONS 

In this paper, we have analyzed the effects of electron- 
electron interactions on the electronic properties of mag- 
netic waveguides formed in suitable inhomogeneous mag- 
netic fields in graphene. When there are two parallel 
lines along which the magnetic field vanishes, a pair of 
counterpropagating snake states can be formed, which 
are ideal unidirectional (chiral) channels similar to the 
edge states in quantum Hall bars. We have studied the 
case of a smooth magnetic field profile across the wire, 
but similar results are expected also for other profiles, 
e.g. for piece-wise constant fields. Employing a combi- 
nation of analytical methods for the band structure and 
for the eigenfunctions (which were checked against ex- 
act diagonalization results) with a thermodynamical ap- 
proach, we have obtained a closed result for the non- 
universal TLL parameter, see Eq. (|32p and Fig. [5l This 
parameter then determines the power-law exponents ap- 
pearing in many observables of interest, in particular in 
the energy-dependence of the tunneling density of states. 
Quite remarkably, we have uncovered that the snake or- 
bits impose a periodic superstructure that breaks Galilei 
invariance of the resulting wire, and modifies the com- 
monly used estimate for the TLL parameter g. The fi- 
nite 7?,-term in Eq. (P5)) refiects this physics. We ex- 
pect this correction to affect also edge states in quantum 
Hall bars. Typical values for g found here are compa- 
rable to the values reported for semiconductor quantum 
wires and single-wall carbon nanotubes. We thus expect 
that non- Fermi liquid behavior should be sufficiently pro- 
nounced to be observable in systems based on inhomoge- 
neous magnetic fields. 

There are two main advantages regarding experimental 
observability in our system when compared to previous 
TLL realizations. First, the TLL parameter g can be 
tuned over a significant region by sweeping the strength 



of the homogeneous part of the magnetic profile. This 
is illustrated in the main panel of Fig. [S] Second, the 
unavoidable presence of disorder is not expected to affect 
the TLL behavior, since right- and left-moving electrons 
are spatially separated. This may allow, for the first time, 
to experimentally study the physics of an ultraclean TLL 
state. 

While there are similarities to the physics of quantum 
Hall edge states, the TLL state discussed here is quite 
different from the chiral Luttinger liquid discussed in the 
context of the fractional quantum Hall effect 4^ The g 
parameter in the latter case is fixed by the bulk filling 
factor, while here it is non-universal and tunable. From 
a conceptual point of view, the quantum Hall situation is 
also more intricate because of the coupling of edge states 
to bulk states4S. Such complications are absent for the 
TLL state discussed in this paper. To conclude, we hope 
that our work motivates experiments in this direction. 
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NUMERICAL DIAGONALIZATION 

In this Appendix, we provide some details on how 
the numerical diagonalization mentioned in Sec. IHBI 
has been implemented. From Eq. ([6]), we first observe 
that exact eigenfunctions 0(.t) of Eq. ([3]) behave as 

~ exp — ^ for a; — > oo. For example, the exact zero- 
energy state to Eq. © is given as (e-?''/'*+''o?'/2+«?, 0)^ 
up to a normalization factor. 

However, it is clear that accurate eigenvalues ej, only 
require to maximize the overlap between the approxi- 



mated wavefunction and 4>{x). For convenience, we thus 
use the complete and orthonormal oscillator basis 

where Hn are Hermite polynomials. Matrix elements 
of the operator H"^ in Eq. ([3]) in this basis read (6o — 

Bq/ujb) 

UJj^^{ipn\H^\ipn') = (2n + 1 + K^)(5„„' 

-(^„|2fc(r - boO - r.H"'' - bo)\Vn') , 

where we use that for integer 7 and even n + ?i' + 7 (F 
denotes the Gamma function) 



ifnlClVn' 



— — p j— H" [n/2] [n72] , , m+m' 

2"+" n \n'\ Y^ ■^ / 1 \ 



m—O m'—O 



p. l+n+«'+7 „^„^/) 
V - 

m\ m'\ [n — 2m)\ {n' — 2to')! 

where {^n\£,^ \y^n') = otherwise. The symbol [n] de- 
notes the largest integer smaller or equal to n. Upon 
carrying out standard diagonalization for < n, n' < 30 
basis functions, we obtain the energy dispersion Sk in 
Figs. [T] and [3] for j/ = 2 and v = 3, respectively. To 
numerical accuracy, all levels are found independent of 
Tz = ±1, although the effective potentials in Eq. ^ dif- 
fer considerably. Only the zero-energy level Ek — in 
Fig. [3] (there is no zero-energy level in Fig. [1]) belongs 
purely to the upper pseudo-spin component t^ = +1 for 
the valley point K chosen here. With the numerical diag- 
onalization, we also obtain the eigenfunctions ipkix) and 
Xk{x) of the Dirac Hamiltonian ^. In Figs. [Handel we 
show the resulting density profiles for v = 2 and i^ = 3, 
respectively. All figures nicely confirm the overall picture 
developed in Sec. IIIBI 
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